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LAUNCH WINDOW ANALYSIS IN A NEW PERSPECTIVE WITH 

EXAMPLES OF DEPARTURES FROM EARTH TO MARS 

By Joseph R. Thibodeau III and Victor R. Bond 
Manned Spacecraft Center 

SUMMARY 


A technique relying upon the effects of gravitational harmonics for parking-orbit 
alinement is analyzed. This technique is applied to the problem of injection from an 
assembly parking orbit around the Earth to an escape trajectory toward Mars. The 
initial orientation of the assembly orbit is chosen so that resulting nodal and apsidal 
perturbations will maintain continually a nearly ideal alinement with the required de- 
parture asymptote during succeeding days of the departure opportunity. Thus, the 
motion of the parking orbit is synchronized with the motion of the escape asymptote. 
Injection velocity penalties caused by misalinement are greatly reduced or eliminated. 

Parking orbits that track the escape asymptote have moderate-to-high apoapsis 
altitudes, typically in a range from 2000 to 8000 nautical miles. This method may be 
useful for planning round trip interplanetary missions. In relation to the high apoapsis 
altitudes, the method could be used for defining parking orbits at a target planet by 
considering the departure window for the return trip to Earth. 


INTRODUCTION 


Conditions for an ideal minimum AV impulse to transfer from a circular parking 
orbit to a hyperbolic escape orbit occur when the orbits are coplanar and the impulse is 
applied at periapsis of the escape hyperbola. An additional condition for an elliptical 
orbit is that the periapsis position vectors of the parking orbit and the escape hyperbola 
must be colinear. For preliminary analysis of an interplanetary mission, it is often 
economical and judicious to assume implicitly that these ideal alinement conditions will 
exist at the time of injection and then to calculate the ideal or characteristic velocity 
of an escape trajectory. 

Unfortunately, these conditions seldom exist when dynamic effects are simulated 
in an analysis of the injection window. First, the oblateness of the planet perturbs the 
ascending node and the argument of periapsis of the orbit. The secular effects of these 
perturbations are that the ascending node and argument of periapsis change linearly 
with time. Second, the hyperbolic excess velocity vector is also changing with time. 
These motions are independent and invariably result in misalinement of the orbital 



plane and the departure trajectory. The escape maneuver no longer can be applied at 
periapsis, and, further, a large plane change may be necessary. Substantial, and 
sometimes severe, velocity penalties may be required to accomplish the required 
turning. 

Many authors have studied this problem (refs. 1 and 2). All of these authors have 
approached the problem by choosing the inclination, eccentricity, and periapsis radius 
of the parking orbit, and then waiting for the ascending node and argument of periapsis 
to change into a configuration such that the velocity impulse required to inject into the 
escape orbit is a minimum. This procedure results in long time periods during a de- 
parture opportunity when injection is impossible because of spacecraft fuel limitations. 

Efficient propulsive techniques have been devised to accommodate this problem 
and to buy more time for scheduling the injection maneuver. Even though these tech- 
niques can reduce the velocity requirements dramatically, there still may be relatively 
long periods during a departure opportunity when injection is impossible because of a 
large plane-change velocity penalty. 

A technique that predicts the parking- orbit inclination and eccentricity required 
to keep the motion of the parking orbit and the escape orbit in near synchronization is 
presented in this paper. The initial orientation of the assembly orbit is chosen so that 
the resulting nodal and apsidal perturbations will shift the original orbit so that it con- 
tinually maintains a nearly ideal alinement with the required departure asymptote dur- 
ing succeeding days of the departure opportunity. Plane-change velocity penalties 
caused by orbit-plane misalinement are reduced greatly or are eliminated. 

Using assembly orbits that track the escape trajectory is an interesting mecha- 
nism for buying more time for scheduling the injection maneuver. Ideally, this ap- 
proach could also afford the most flexible schedule for launch and orbital assembly of 
the spacecraft and the propulsion systems. Unexpected difficulties can delay the launch 
of a vehicle or the assembly of several vehicles in Earth orbit. In addition to built-in 
holds in the assembly schedule, delays can be accommodated by letting the injection 
date slip to a later time within the injection window. Slipping the injection date no 
longer has a consequence because there is an opportunity to initiate the escape maneu- 
ver during every revolution of the parking orbit. Therefore, in the ideal situation, this 
technique offers flexibility as a principal advantage. 

Booster performance and the geographic location of the launch site largely dictate 
the characteristics of the assembly orbit; the general problem is to maximize the pay- 
load that can be delivered to the assembly orbit. Ordinarily, the spacecraft and pro- 
pulsion systems would be launched directly into an assembly orbit, which not only 
maximizes the payload but also affords a timely schedule for launching several vehicles. 
These considerations lead to the principal (but unproven) disadvantage of the technique 
presented in this report. That is, the characteristics of assembly orbits that track the 
escape asymptote may or may not be compatible with the characteristics dictated by 
booster performance. This problem cannot be answered until the characteristics of 
orbits that are dictated by the two methods are discovered. The ideal situation would 
be that the characteristics of the assembly orbit dictated by the two methods would be 
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the same or at least compatible. This report investigates the following two ques- 
tions and presents a partial resolution of the problem. 

1. What are the orbital elements of the parking orbits that track the escape 
asymptote and permit a departure at any time during the launch opportunity? 

2. How long does the motion of the orbit remain synchronous with the motion of 
the departure asymptote? 

The authors thank Ellis W. Henry and Gregory A. Zambo, NASA Manned Space- 
craft Center, Houston, Texas, for contributing their work on the multiple-impulse 
heliocentric transfer trajectories. 


SYMBOLS 


ALFA 
APO ALT 

a 

C 

CIR DV 

c 


right ascension of excess hyperbolic velocity vector, deg 
apoapsis altitude, n. mi. 
semimajor axis 

constant defined by equation (22) 
circular-to-elliptical transfer AV, fps 
auxiliary constant for the differential 


DELTA 

ECC 

EMOS 

e 

f 

g 

h 

INCL 

i 


declination of excess hyperbolic velocity vector, deg 
eccentricity (FORTRAN symbol in tables in, IV, and V) 
Earth mean orbital speed 
eccentricity 

auxiliary variable for the differential 
auxiliary variable for the differential 
auxiliary variable in equation (16) 

inclination, deg (FORTRAN symbol in tables III, IV, and V) 
inclination, deg 

critical inclination equal to 63. 43 deg 
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k 

NODE 

n 

OMEGA 
PC ANG 
PERIOD 

R B 


oblateness coefficient 

indexing parameter (k = ± 1) 

right ascension of ascending node, deg 

mean motion 

argument of periapsis, deg 
plane-change angle, deg 
orbital period, hr 
equatorial planet radius 


r 

IT 


radius of periapsis 


SMA 
SUM DV 
TEI DV 
TJD 
TMI DV 
TRUAN 
t 

Vco 

VMAG 
V PERI 
Z 


semimajor axis, ft 

sum of CIR DV and TMI DV 

injection AV at the time of ideal alinement 

elliptical orbit departure Julian date, days 

trans-Mars injection AV, fps 

true anomaly of departure maneuver 

time 

excess hyperbolic velocity 

excess hyperbolic velocity magnitude, fps 
velocity at periapsis, fps 
variable defined by equation (20) 

variables introduced for solution of cubic equation (A2) 
right ascension of the V^ vector 


cKq, aj, a g coefficients of £; in cubic equation (A2) 
AV impulsive-velocity increment 
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6 declination of the V OT vector 

r] true anomaly of hyperbolic escape asymptote 

p. gravitational parameter 

2 

£ trigonometric variable equal to tan i 

cr auxiliary variable used in defining £2^ 

£2^ ascending node as computed from escape asymptote 

£2 secular rate of ascending node computed from planetary oblateness 

s 

£2^ secular rate of ascending node computed from motion of escape asymptote 

argument of periapsis as computed from escape asymptote 

ci> secular rate of argument of periapsis computed from planetary oblateness 

s 

w M secular rate of argument of periapsis as computed from motion of escape 

asymptote 

ANALYSIS 


The method that will be presented is similar to the one used for alining a parking 
orbit about an oblate planet with arrival and departure hyperbolic excess velocity vectors 
as presented in reference 3. In that reference, alining an orbit at an oblate planet re- 
quires that the total change in the ascending node and argument of periapsis of the park- 
ing orbit be equal to the corresponding quantities of the escape orbit. The method used 
here requires that the ascending node rate and argument of periapsis rate of the parking 
orbit match instantaneously with the corresponding rates of the escape orbit. 


Derivation of Synchronous i and e 

The requirement that the two orbits be synchronized at any time t may be ex- 
pressed mathematically as 


« s (t) = t) 

(i) 

« s (t) = vjt) 

(2) 
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It will be assumed in this derivation that the radius of periapsis r^ is specified 
and that the hyperbolic excess-velocity vector (V B , a^, 6^) and its rate (V^, a^, 
6) as functions of time are known. 

OO' 

From reference 4, J2 and o> are found from 
’ s s 


-3nJ„R * 

fi = ^ (cos i) 


2a‘ 


! (l-e 2 )‘ 


0) 


- 3nJ 2 R B 


2a 


J (l-e 2 )' 


(§ sin2i - 2 ) 


(3) 


(4) 


It must be noted that these are secular perturbations (as implied by the subscript s). 

Because i, a, and e are not secularly perturbed, Cl and u> g are constants, when 
i, a, and e are known. ° ° 

From reference 5, the ascending node and argument of periapsis of the hyperbola 
are given at any time by 


fi OO = Q 00 + ” + 


or 


Cl 


OO 


a 


a 


where 


a - 


. -1 
sin 


' tan O 

i tan i / 


(5) 


( 6 ) 


(7) 


As shown in reference 5, equation (5) is used for a maximum periapsis declination 
and equation (6) for a minimum periapsis declination. 
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The argument of periapsis of the hyperbola is given by 


w = -77 + tan 


-1 


tan 6„ 


sin i cos^ - a ^ 


( 8 ) 


where 


7 ] = cos 

The time derivatives of £2^ and are found by differentiating equations (5) and (6) 
and equation (8). These results are 


- 1 + 



0) 


S2 = 

OO 


<*oo + f5 oo 


( 10 ) 


^co = CV oo + g 6 o 


where 


f = 


cos cos a tan i 

CO 


( 11 ) 


(12) 


g = 


-k 

sin i cos a 


(13) 


c 


o 2 V r 
2 COS 77 °° 77 

sin 77 p. 


(14) 


and k = ± 1, depending upon whether maximum or minimum periapsis declination is 
desired. 
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From equations (1) and (2) 


fi & 
oo s 


U> CO 
oo g 


By using equations (3), (4), (10), and (11) in equation (15) 


a + f6 

00 00 , COS 1 

= h = 


,, t 5 . 2. _ 

cV co + e 6 oo 7 T sin i - 2 


(15) 


(16) 


Note that equation (16) is independent of the eccentricity e. Thus, it may be solved by 
iteration for the inclination i, because and are known. An interesting par- 
ticular solution of equation (16), when V M = 0, is given in the appendix. This solution 

is analytic and is useful in providing initial guesses for the numerical solution of 
equation (16). 

Three restrictions on the inclination should be stated. The first restriction is 

that 


it - Ul > i > I sl (17) 

which is seen by inspection of equation (7). The second restriction is that if 


O < 0 

oo 



OO 


then 

. ^ 7 J 
X< 2 


and if 

Ro> 0 


then 


(18) 


which is seen by inspection of equation (3). The third restriction is that if 


u> < 0 

OO 

then i L < i < 77 - i L 
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r 


and if > 0 

then 0 < i < (19) 


where i L 



63.43°, 


which is seen by inspection of equation (4). 


Violation of the restriction of equation (17) will cause the iterative solution of 
equation (16) for i to fail immediately. The restrictions of equations (18) and (19) 
must be tested before the determination of the eccentricity is begun. 


When the inclination has been found, the eccentricity may be found from equa- 
tions (1) and (3). 


S2 




3n J 2 R B 


cos 1 cos 1 


2a‘ 


! (l - e 2 )‘ 


= z 


(20) 


Because a = r /(I - e) and 

7T ' 



equation (20) becomes 


Z _ (1 - e) 
C 


3/2 


(1 +e) £ 


( 21 ) 


where 


C = - 


R 2 
2 J 2 K B 


2 nr 

Jr 7 

1 7T 


( 22 ) 


Equation (21) may now be solved by iteration for the eccentricity e and, in fact, is 
identical to the iteration performed in reference 3. Note from equation (21) that the 
quantity Z/C must always be positive for a solution to exist. The requirement that 
Z/C be positive is simply a restatement of the restriction given by equation (18). 
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The computation of the inclination and eccentricity can now be summarized. 

1. Given V M , at time t, r , and Z, solve equations (16) for the 

inclination. 

2. Compute £2^ and from equations (10) and (11); then test this inclination 
to see that it meets the restrictions of equations (18) and (19). 

3. Now solve equation (21) for the eccentricity. 


PROCEDURE 


The techniques discussed in this report have been programed in FORTRAN V for 
the UNIVACT108 digital computer. Two approaches were investigated. In the first 
approach, instantaneous choices of the orbital elements were calculated as functions 
of time. By making the proper choices of inclination and eccentricity (ideally, updating 
these elements continuously by small midcourse maneuvers), the orbit would be forced 
to maintain ideal alinement. Although the individual midcourse maneuvers are small 
(from 80 to 100 fps to as large as 500 fps every 5 days), their sum rapidly accumulates; 
consequently, this approach was abandoned. 

In the second approach, the parking orbit is not forced to remain ideally alined. 
Instead, the orbit is required to be ideally alined at some time during the injection 
window and to be approximately alined during the remainder of the injection window. 

For this approach, there are no guidance and targeting corrections; it is assumed that 
all the accumulated errors are removed at the time of the injection maneuver (or ma- 
neuvers as in the case of a multiple-impulse departure sequence). Thus, the orbit 
merely drifts and, in a general sense, follows the motion of the escape asymptote. 

The orbit is no longer ideally alined. The injection maneuver may again include an 
off-periapsis thrust as well as a plane change. The following are important differences. 

1. The orbit can be alined during any portion of the injection window. 

2. The plane-change misalinement is greatly reduced or eliminated during the 
remainder of the window. 

The second approach has been used to determine the effective Earth-orbital depar- 
ture windows for two types of Mars round trip stopover missions: the 1981 Mars stop- 
over with inbound Venus swingby and the 1986 direct Mars stopover with near-Hohmann 
transfers. The heliocentric profiles for these missions are shown in figures 1 and 2. 

The Mars mission flight plans were calculated by using a patched-conic trajectory 
model. Successive conic trajectories were pieced together to form a reference mission 
or flight plan; a linear weighted least- squares iteration technique was used to produce 
the launch dates, flight ti-mes, and impulsive characteristic velocity requirements for 
flight plans that exhibit a minimum in the sum of the velocity impulses. Patched-conic 
solutions were calculated at 5-day steps through 100-day Earth departure periods in 
1981 for the swingby mission and in 1986 for the direct mission. The hyperbolic excess 
velocity vectors at Earth departure were then used to compute the characteristics of 
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Departure Hyperbola 




(?) Mars departure 
(?) Earth arrival 

Figure 2. - Illustration of the conjunction 
class mission and the parking orbit 
about Mars. 


Figure 1.- Illustration of the 1981 Mars 
stopover /Venus swingby mission. 


the parking orbits at Earth and the 
injection-velocity requirements for the 
Earth escape maneuver. 

The Earth escape maneuver is mod- 
eled by a single, minimum AV, impulsive 
thrust that includes a plane change as well 
as an off-periapsis thrust component, as 
shown in figure 3. The general formulation 
for the injection AV is given in refer- 
ence 5. The AV is minimized by a one- 
dimensional search on the parking orbit 
true anomaly. 


Norpial 
to orbit 
plane 



Figure 3.- Geometry of single- impulse 
injection maneuver. 


RESULTS 


The motion of the escape asymptote during the course of a launch opportunity 
dramatically affects the solutions that can be obtained for the tracking orbit. There 
is no question that solutions can be found. Whether these solutions are useful or desir- 
able is the question that must be investigated. 

A desirable solution may have restrictions on the range of acceptable orbital 
inclinations and apoapsis altitudes. Because these parameters are the dependent vari- 
ables of a deterministic system of equations, there is no discrete control over them. 
Two methods to control the inclination and eccentricity of the tracking orbit can be 
investigated. 

1. Trajectory shaping can be used to control the orientation and motion of the 
escape asymptote. This approach is demonstrated by the use of two- and three- impulse 
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2. Given the time history of the escape asymptote during the course of a launch 
opportunity, solutions can be obtained at discrete times during the launch opportunity. 
This approach amounts to shifting the epoch of the solution (e. g. , the time of ideal 
alinement) and provides a mechanism for "fine tuning" the orbital elements. 

Both approaches have been applied in an analysis of the 1981 Venus swingby and 
the 1986 direct minimum- energy Mars stopover missions. The results for the two- 
impulse trajectories are discussed in the following paragraphs, and the results for the 
three- impulse trajectories are discussed in a later section. 

The analysis for the 1981 Venus swingby mission is included as an example to 
show the effect of shifting the time of ideal alinement (e. g. , the epoch of the solution). 
The results presented for this mission are typical. The orbits display high apoapsis 
altitudes and both posigrade and retrograde inclinations. The solutions for the orbit 
inclinations and apoapsis altitudes are found to be slightly different when the solutions 
are evaluated at different epochs. Typically, three solutions are found at any epoch. 
The reasons for this are discussed in the appendix. The orbits track for approximately 
50 percent of the total available window. 


The 1981 Mars Stopover With Venus Swingby 

The 1981 Mars stopover with a Venus swingby mission consists of a direct inter- 
planetary transfer from Earth to Mars, a short orbital stopover at Mars, and a return 
transfer from Mars to Earth that is also a free flyby of Venus. The direct transfer 
from Earth to Mars can be either a type I (heliocentric transfer angle <180°) or a 
type II (heliocentric transfer angle > 180°) trajectory. Both transfer trajectories have 
the same arrival dates at Mars. The Earth-departure dates for these trajectories are 
separated by 40 to 60 days. The result is that there are two separate Earth -departure 
opportunities, one in November of 1981 and the other in January 1982. A summary of 
the characteristics of these missions is presented in tables I and II. The stay time at 
Mars for all swingby missions analyzed in this report is 60 days. 

The Earth -departure velocity contour for this mission displays a characteristic 
double minimum, and there is a rise in the departure velocity at the transition from 
the type I to the type n solutions. The double minimum in the Earth- departure velocity 
contour prevents the use of the total available launch window. Tracking orbits can 
follow the motion of the escape asymptote only if the motion is continuous. Thus, the 
type I and type n departure windows must be considered separately. (Using broken- 
plane heliocentric trajectories with an intermediate impulse applied between Earth and 
Mars can eliminate the discontinuity. This will be discussed later. ) 

The Earth-departure velocity requirements for the type II window are shown in 
figures 4(a), 4(b), and 4(c). The effect of shifting the epoch of the solution is found by 
comparing these figures. The time of ideal alinement is indicated in the figure where 
the curves converge and is marked by asterisks and a crosshatched region at the bottom 
of the grid. Three solutions are plotted in figure 4(a). These solutions have inclina- 
tions and apoapsis altitudes of 36° and 29 300 nautical miles, 99° and 5800 nautical 
miles, and 104° and 5800 nautical miles. 
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39*, 17 000 n. mi. 


16 x 10 3 



(a) Aline ment date of October 31, 1981 
(Julian date 2 444 908, 65). 


16 x io 3 



o 1 > ' ‘ 1 J 

4890 4900 4910 4920 4930 4940 4950 

Date of Earth departure, days after Julian date 2 440 000 


(b) Alinement date of November 5, 1981 
(Julian date 2 444 913. 65). 



(c) Alinement date of November 10, 1981 
(Julian date 2 444 918. 65).. 

Figure 4. - Earth- departure velocity 
requirements for the 1981 Venus 
swingby with Earth-to-Mars trans- 
fer angles greater than 180°. 


The fourth and bottom curve displays 
the ideal or characteristic velocity for de- 
parture from a 262-nautical-mile altitude 
circular-parking orbit. This curve is 
simply the scalar difference between the 
periapsis speed of the escape hyperbola and 
the circular speed of the 262-nautical-mile 
altitude parking orbit, and it represents 
the absolute minimum injection AV. The 
bottom curve thus provides a standard of 
comparison to measure the performance of 
different solutions. The bottom curve, used 
as a basis for comparison, is indicative that 
the 104° orbit, which appears to maintain 
alinement best, tracks almost ideally the 
first 15 days before gradually drifting out of 
alinement. Likewise, the 36° inclined orbit 
appears to track the first 15 days but rapidly 
diverges during the remainder of the launch 
window. 


If the time of alinement is shifted by 
10 days, there are two solutions as shown 
in figure 4(b); one with an inclination of 97° and apoapsis altitude of 5500 nautical miles, 
and the other with an inclination of 39° and 17 000- nautical-mile altitude. 


The apoapsis altitude of 17 000 nautical miles for the 39° solution is considerably 
lower than the corresponding solution of the 29 000 -nautical -mile altitude and 36° incli- 
nation, as shown in figure 4(a). 


If the time of alinement is again shifted by 10 days, two more solutions, which 
are shown in figure 4(c), are found. The solutions are at 41° inclination with an 
11 000-nautical-mile apoapsis altitude and at 127° inclination with a 7 3 00- nautical- 
mile altitude. When the epoch is shifted forward through the launch period, the 
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apoapsis altitude decreases, and different solutions are found at each new epoch. The 
curves of figure 4 are typical of solutions that can be found. If an artificial limit is 
drawn to represent a maximum permissible injection AV of 12 600 fps (represented 
by the broken line in fig. 4(c)), it can be seen that there are 30 days during the type II 
window in which the ideal injection AV is less than 12 600 fps. The solution with the 
41° inclination provides 22 days during which the injection AV is less than 12 600 fps. 
This solution tracks approximately 67 percent of the available launch window. 

The time histories of the solutions in figure 4 are presented in tables 111(a), m(b), 
and in(c). The plane change and off-periapsis thrust requirements, the orbital ele- 
ments, and other parameters that describe the solution are presented. These param- 
eters are included for future reference and possible comparison with other methods. 

In table HI, the time of ideal alinement is indicated when the off-periapsis thrust 
angle (TRUAN) and the plane-change angle (PC ANG) are simultaneously zero. Two 
values of AV (TMI DV and SUM DV) are shown in table in. The TMI DV is the vec- 
tor AV required for a single-impulse injection maneuver to escape from the elliptical 
tracking orbit into a trajectory toward Mars. The SUM DV considers the additional AV 
that is required to transfer from the initial 262-nautical-mile altitude circular parking- 
orbit to the elliptical orbit. The circular-to- elliptical orbit transfer AV is printed 
under the heading CIR DV. 

The time history of the motion of the escape asymptote is also shown in table in. 
The excess hyperbolic velocity vector is printed in terms of speed in feet per second 
(VMAG), and right ascension (ALFA) and declination (DELTA) in degrees. The decli- 
nation is measured with respect to the 
equatorial plane of the Earth. The right 
ascension is measured from an inertial 
X-axis through Aries; this is also the 
system in which orbital elements are 
measured. 

The Earth-departure velocity require- 
ments for the type I launch window are 
shown in figures 5(a), 5(b), and 5(c). The 
time histories of the solutions are pre- 
sented in tables IV(a), IV(b), and IV(c). 

Three solutions are shown in figure 5(a); 
one is posigrade with a 46° inclination and 
4500-nautical-mile apoapsis altitude, one 
is nearly polar with an 81° inclination and 
4 8 63 -nautical-mile altitude, and one is 
retrograde with a 147° inclination and 
5815-nautical-mile altitude. An examina- 
tion of figures 5(b) and 5(c) indicates the 
effect of shifting the epoch forward through 
this launch window. The inclinations re- 
main nearly the same, and the altitudes of 
the solutions increase. 


16 x 10 3 



(a) Alinement date of December 20, 1981 
(Julian date 2 444 959. 31). 

Figure 5.- Earth- departure velocity 
requirements for the 1981 Venus 
swingby with Earth-to-Mars trans- 
fer angles less than 180°. 
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(b) Alinement date of December 25, 1981 (c) Alinement date of December 30, 1981 

(Julian date 2 444 964. 31). (Julian date 2 444 969. 31). 

Figure 5.- Concluded. 


Also in figure 5(a), an artificial limit can be drawn to represent a maximum 
permissible injection AV of 12 600 fps (represented by the broken line in fig. 5(a)). 
There are 35 days in the type I window during which the ideal injection AV is less than 
12 600 fps. The solution with the 46° inclination provides 30 days during which the ve- 
locity requirements are less than 12 600 fps. The 46° inclined orbit permits departure 
any time (i. e. , at each periapsis passage) during 86 percent of the available launch 
window. 

The injection maneuver is a sequence of two maneuvers, the CIR DV and the 
TEI DV. Often these are of nearly equal magnitude. There are several methods of han- 
dling the implied logistical problems of launch, Earth-orbit assembly, and scheduling of 
these two maneuvers. For example, if the altitude of the solution is low enough, the 
spacecraft could be launched directly into a tracking orbit; if the altitude is high, the 
spacecraft could be launched into an intermediate phasing orbit and the CIR DV could be 
scheduled after the launch and assembly and before the opening of the injection window. 

The CIR DV is introduced to form a relative basis for comparison of different 
solutions and it may or may not be scheduled discretely in an operational sequence. 


The 1986 Direct Minimum-Energy Mars Stopover 

The direct minimum- energy Mars stopover mission consists of a near-Hohmann 
transfer from Earth to Mars, a stopover at Mars of 300 to 500 days, and a near- 
Hohmann return transfer from Mars to Earth. The total trip times for these missions 
are approximately 750 to 1000 days. 

The Earth-departure opportunity in 1986 allows 110 days during which the ideal 
injection AV (for departure from a 262-nautical-mile altitude circular orbit) is less 
than 12 600 fps. Both type I and type II transfers simultaneously take on minimum 
injection AV’s during this departure period, and trajectories with 180° Earth-Mars 
transfer angles do not present a problem as in the 1981 Venus swingby mission. 
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The 1986 direct mission opportunity contains unique problems. The type I trans- 
fers (transfer angles less than 180°) have rather high departure asymptote declinations. 
The asymptote declination is greater than 40° for more than half the departure oppor- 
tunity and is as high as 62° in the early part of the opportunity. It is not possible to use 
an assembly-parking orbit with an inclination less than these declinations without plan- 
ning for significant plane- change velocity penalties. A parking orbit that will track 
these asymptotes without a plane-change requirement must have a high orbital 
inclination. 

The use of type II trajectories eliminates the problem of high asymptote declina- 
tions. In fact, for most of the type n window, the asymptote declination is within 5° of 
the Earth's equatorial plane. However, the type II trajectories introduce a new prob- 
lem. The motion of the escape asymptote is such that only parking orbits with high 
retrograde inclinations will track them for very long. 

The solutions for the type I window are shown in figure 6(a) and table V(a). Two 
solutions are shown in figure 6(a), one at 78° inclination with a 5 3 00- nautical-mile 
apoapsis altitude and the other at 55° inclination with a 1200-nautical-mile apoapsis 
altitude. The 78° inclined orbit provides 50 days during which the velocity require- 
ment is less than 12 600 fps. The 55° inclined orbit provides almost 30 days during 
which the departure velocity is less than 12 600 fps. As can be seen in figure 6(a), the 
orbits track moderately well for one-third to one-half the total available injection win- 
dow. These two solutions are shown in more detail in table V(a). 

The solutions for the type II window are shown in figure 6(b) and table V(b). 

Three solutions are shown in figure 6(b), and they all display very high apoapsis 




(a) Type I heliocentric transfers. (b) Type n heliocentric transfers. 

Figure 6. - Earth- departure velocity requirements for the two-impulse direct 

mission in 1986. 
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altitudes: 33° and 21 600 nautical miles, 97° and 9350 nautical miles, and 126° and 
13 000 nautical miles. The 33° inclined solution allows a 30-day departure period dur- 
ing which the departure velocity is less than 12 600 fps. The 97° inclined solutions 
yield 40 days in which the departure velocity is less than 12 600 fps. These two solu- 
tions track for about 30 to 40 percent of the total available injection window. The third 
solution at 126° provides nearly 60 days during which the departure velocity is less 
than 12 600 fps, or about 75 percent of the total available injection window. These solu- 
tions are shown in more detail in table V(b). 


Multiple- Impulse Heliocentric Transfer Trajectories 

The Mars orbit, which is slightly elliptical and inclined to the Earth’s orbital 
plane, precludes the possibility of a true Hohmann trajectory. More important, the 
injection geometry and AV requirements change during the course of the departure 
opportunity. There is flexibility for controlling these requirements by use of multiple- 
impulse heliocentric trajectories, particularly three-impulse trajectories in the tran- 
sition from type I to type II solutions. 

For single-impulse ballistic trajectories, the heliocentric paths are represented 
by conic arcs that are solutions of Lambert’s problem. The conic arc will have the 
Sim at its focus and will lie in a plane containing the Sun and the position vectors of the 
Earth at the time of departure, and the target-planet Mars at the time of arrival. When 
the difference in heliocentric longitude between the Earth at departure and the target 
planet at arrival approaches 180°, the plane containing the heliocentric transfer be- 
comes highly inclined to the ecliptic and, in the limit, becomes polar. Although the 
inertial velocity has not increased for these limiting solutions, the relative Earth de- 
parture AV increases tremendously (on the order of one Earth mean orbital speed 
which is equal to approximately 100 000 fps). The effect of the Earth's orbital speed 
around the Sun is lost and nearly all the inertial velocity must be obtained through the 
spacecraft propulsion system. 

High-Earth-departure AV’s can be avoided by permitting an intermediate impulse 
on the Earth-to-Mars trajectory. With the intermediate impulse, the heliocentric tra- 
jectory is broken into two legs. The legs of the new trajectory intersect and have a 
common position and time in space. In the general problem, four independent variables 
are introduced: three coordinates to describe the position of the impulse and one to 
describe the time of the impulse. The problem ceases to be deterministic, and many 
solutions can be found. 

Three kinds of three-impulse trajectories were investigated, one for the 1981 
Venus swingby and two for the 1986 direct mission. A minimum-sum AV three- 
impulse trajectory was used for the Venus swingby with no control over the geometry of 
either leg. The method used is presented in reference 6. 

A minimum-sum AV three-impulse trajectory was also used for the 1986 direct 
mission. For these trajectories, the geometry of Earth-departure and the Mars- 
arrival trajectories was controlled. In one kind, the Earth- departure trajectory is 
constrained to lie in the ecliptic plane; in the other, the Mars-arrival trajectory is con- 
strained to lie in the Mars -orbital plane. 
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The 1981 Venus Swingby With Three-Impulse Outbound Trajectory 


The 1981 Venus swingby mission is somewhat improved using three-impulse tra- 
jectories. The best orbit, shown in figure 7 and table VI, tracks for nearly 45 days 
with the AV for Earth departure less than 12 600 fps or for more than half the total . 
available injection window of 70 days. This orbit has an inclination of 75° and an 
apoapsis altitude of 1900 nautical miles. Other solutions found for this mission are 
summarized in table VII. 



Figure 7.- Earth- departure velocity re- 
quirements for the 1981 Venus swing- 
by with minimum 2AV three-impulse 
trans-Mars trajectories. 


One may note that there still is a dis- 
continuity in the velocity requirements for 
Earth departure and that the orbit does not 
track the full injection window in figure 7 
partly because of the discontinuity. It is 
not understood why this discontinuity exists, 
because the Hohmann spike has been elim- 
inated. The clue to the problem is that, in 
the mission analysis program, the Mars 
orbit stay time was forced to be 60 days and 
the Venus flyby was forced to be a free flyby 
with a positive flyby altitude. The result is 
that the return leg from Mars to Earth is 
nearly always the same, and changing the 
outbound flight time from Earth to Mars 
accounts for the variability in the Earth de- 
parture date. (This is a result and not a 
restriction in the model because the flight 
times were independent variables chosen to 
minimize the sum of all the velocity im- 
pulses. ) A similar result can be seen in 
the two- impulse trajectories shown in 
tables I and II. These missions are mini- 
mum sum AV mission plans with fixed 
Mars stay times. Note that the Venus swing- 
by date is always the same. 


The 1986 Direct Mission With Three- Impulse Trajectories 

The results for the 1986 direct mission using three-impulse trajectories are 
shown in figures 8 and 9. Constraining the Earth -departure asymptote to the ecliptic 
plane for the intermediate impulse trajectory in 1986 has the effect of arresting the 
motion of the escape asymptote. The end effect is that very slow regression rates are 
needed to track the slow- moving asymptote. The result is that extremely high apoapsis 
altitudes are found, as is shown in figure 8. One orbit at 18° inclination tracks ex- 
tremely well, tracking nearly 90 days with the departure velocity below 12 600 fps. 

The high apoapsis altitudes are similar to those in the type n trajectories shown in 
figure 6(b). 
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Earth departure velocity, fps Earth departure velocity, fps Earth departure velocity, (ps 


I 



6480 6500 6520 6560 6560 6580 6600 


Date of Earth departure, days after Julian date 2 440 000 

(a) Date of alinement at 2 446 510. 


L5 x 10 3 



Date of Earth departure, days after Julian date 2 440 000 

(b) Date of alinement at 2 446 520. 
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6480 6500 6520 6540 6560 6580 6600 

Date of Earth departure, days alter Julian date 2 440 000 

(c) Date of alinement at 2 446 530. 

Figure 8. - Earth- departure velocity- 
requirements for the 1986 direct 
Mars mission using three-impulse 
trans-Mars trajectories with the 
Earth- departure asymptotes in the 
ecliptic plane. 



Date of Earth departure, days after Julian date 2 440 000 

(a) Date of alinement at 2 446 545. 



6430 6500 6520 6540 6560 6580 6600 


Date of Earth departure, days after Julian date 2 440 000 

(b) Date of alinement at 2 446 555. 



6480 6500 6520 6540 6560 6580 6600 

Date of Earth departure, days after Julian date 2 440 000 


(c) Date of alinement at 2 446 570. 

Figure 9.- Earth- departure velocity 
requirements for the 1986 direct 
Mars mission using three-impulse 
trans-Mars trajectories with the 
Mars-arrival asymptotes in Mars- 
orbital plane. 
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20 X 10 3 



(d) Date of alinement at 2 446 520. 


Figure 9. - Concluded. 


The results of constraining the Mars- 
arrival asymptote in the Mars- orbit plane 
are shown in figure 9. These results are 
essentially the same as those found in the 
type I launch window shown in figure 6(a). 
There is a somewhat wider variation in the 
choice of orbital inclinations than there is 
for the type I window. The range of apoap- 
sis altitudes is slightly higher than the 
type I window. A retrograde orbit, as 
shown in figure 9(d), of 137° and 
10 400-nautical-mile altitude tracks nearly 
90 days with the departure velocity less 
than 12 600 fps. 

The three-impulse trajectories give 
generally the same result as the two- 
impulse trajectories. The orbits appear 
to track for longer periods of time when 
the three-impulse trajectories are used. 


CONCLUDING REMARKS 


A preliminary injection window analysis was performed for two round trip Mars 
stopover missions: the 1981 Mars stopover with inbound Venus swingby and the 1986 
direct minimum- energy Mars stopover. 

The results showed that by properly choosing the orbital elements of the Earth 
launch and assembly parking orbit, orbits can be found to follow the motion of the re- 
quired trans-Mars departure asymptote. 

For the direct mission using near-Hohmann transfers, the Earth-departure op- 
portunity in 1986 provides nearly 100 days during which the ideal scalar or character- 
istic velocity requirements are less than 12 600 fps for the trans-Mars injection 
maneuver. When parking orbits were alined properly for this mission, they were found 
to track from 50 days to nearly the full 100 days with little or no velocity penalties caused 
by misalinement of the parking orbit and the departure asymptote. The vector-velocity 
requirements can be maintained within one-half of 1 percent of the ideal characteristic 
velocity for periods of as much as 90 days within the departure opportunity. Further, 
solutions are found with orbital inclinations of approximately 30°, 40°, 50°, and 70°. 
These solutions have high apoapsis altitudes, often from 5000 to 20 000 nautical miles. 

For the Mars stopover with the inbound Venus swingby, the Earth -departure op- 
portunity in 1981 provides 70 days during which the ideal scalar AV requirements are 
less than 12 600 fps for the trans-Mars injection. 
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Two solutions for a moderately high elliptical parking orbit maintained alinement 
for periods as long as 40 days during this injection opportunity. The solutions have 
inc lina tions of 53° and 74° with apoapsis altitudes near 2400 nautical miles. In rela- 
tion to the high apoapsis altitudes, the technique may be useful for defining parking 
orbits at a target planet by considering the departure window for the return trip to 
Earth. 


Manned Spacecraft Center 

National Aeronautics and Space Administration 
Houston, Texas, September 27, 1971 
981-30-00-00-72 
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TABLE I.- LAUNCH DATES AND FLIGHT TIMES FOR MINIMUM-SUM AV MARS STOPOVER/VENUS SWINGBY MISSIONS 

LAUNCHED IN 1981 and 1982 


Calendar date 
of launch, 
month/day/yr 

Julian date 
of launch 
(add 2 440 000) 

Outbound 

flight 

time, 

days 

Mars arrival 
date 

(add 2 440 000) 

Orbital 

stay 

time, 

days 

Intermediate 
flight time, 
days 

Venus flyby a 
Julian date 
(add 2 440 000) 

Return 

flight 

time, 

days 

Total 

trip 

time, 

days 

Outbound 

trajectory 

type 

11/19/81 

<-mm 

297 

5225 

20 

144 

5389 

157 

618 

n 

1/10/82 

■H 

233 

5213 

20 

154 

5387 

155 

562 

i 

11/12/81 

4920 

297 

5217 

30 

142 

5389 

157 

626 

- ' * • 

1/8/82 

4978 

228 

5209 

30 

149 

5388 

156 

563 


11/5/81 

4914 

294 

5208 

40 

141 

5389 

157 

632 


1/12/82 

4981 

224 

5205 

40 

143 

5388 

156 

563 


1 1/1 0/81 

4919 

278 

5197 

60 

133 

5390 

152 

623 

n 

12/30/81 

4969 

224 

5193 

60 

137 

5390 

160 

581 

i 

11/9/81 

4918 

267 

5185 

80 

126 

5391 

172 

645 

n 

12/23/81 

4961 

219 

5180 

80 

131 

5391 

165 

595 

i 


a The Venus flyby altitudes vary between 2500 and 2600 n. mi. above the planet’s surface. 


TABLE n. - VELOCITY REQUIREMENTS FOR MINIMUM-SUM AV MARS STOPOVER/VENUS SWINGBY MISSIONS 

LAUNCHED IN 1981 and 1982 


Calendar date 
of launch, 
month/day/yr 

Julian date 
of launch 
(add 2 440 000) 

Orbital 

stay 

time, 

days 

AV of Earth 
departure, a 
fps 

AV of Mars 

arrival, ^ 
fps 

AV of Mars 

departure, ** 
fps 

AV of Earth 

entry, 

fps 

AV 

total, d 
fps 

AV e 
spacecraft, 

fps 

Outbound 

trajectory 

type 

11/19/81 

4928 

20 

11 900 

4100 

11 300 

40 000 

27 300 

15 400 

n 

1/10/82 

4980 

20 

12 800 

5500 

11 400 

40 200 

29 700 

16 900 

i 

11/12/81 

4920 

30 

12 000 

4200 

11 400 

40 000 

27 600 

15 600 

1 B 

1/11/82 

4978 

30 

12 600 

5600 

11 300 

40 000 

29 500 

16 900 

i 

11/5/81 

4914 

40 

12 300 

4500 

11 400 

40 000 

28 200 

15 900 

n 

1/12/82 

4981 

40 

12 400 

5800 

11 400 

40 400 

29 600 

17 200 

i 

11/10/81 

4919 

60 

12 200 

5000 

11 800 

40 300 

29 000 

16 800 

n 

12/30/81 

4969 

60 

12 000 

6390 

11 560 

39 800 

29 950 

17 950 

i 

11/9/81 

4918 

80 

12 300 

5800 

12 600 

39 800 

30 700 

18 400 

n 

12/23/81 

4961 

80 

11 800 

7555 

12 000 

39 700 

31 200 

19 555 

i 


a Parking orbit = 262-n. mi. altitude circular. 
^Parking orbit - 200 by 10 000 n. mi. 
c 400 000-ft entry altitude. 
d Sum of columns 4, 5, and 6. 
e Sum of columns 5 and 6. 
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TABLE m. - THE 1981 VENUS SWINGBY WITH EARTH-TO-MARS TRANSFER ANGLES GREATER 

THAN 180° 

(a) Date of ideal alinement is October 31, 1981 (Julian date 2 444 908. 65) 
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24449)8.65 

4.92 

18-83 

-.bo 

7 26 J • 36 

12141.80 

11 34 1 .Uu 

- 1 7 2 . da 

1 b • 39 

.68 

244492J.65 

6 * 26 

1 b • 3n 

“1*00 

72 f ,4 .46 

l 20 8 2 - 9 u 

1 U 8 j 6 . u u 

-l7l*/8 

1 8 *bb 

2*10 

2444928.6b 

7.60 

11-77 

• b 0 

74 , 1 .b 7 

1 229j . 0 1 

l U442 « Uu 

- 1 7u *89 

22 • 6 1 

3* 9b 

2444933.6b 

8 .9b 

8.2m 

2 -bu 

7 9 2 S • 2 2 

1 2 6 U 3 • 6 6 

1 J 2 7 6 • u J 

- 1 7 u • 1 6 

2 7-66 

b .82 

2444930.65 

10*29 

4 . 7 ? 

3 • (J U 

86 1 9 • 9q 

13490.34 

i U 4 u b • U u 

-lbV.Jb 

33-« J 

7 • b 1 
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TABLE IB. - THE 1981 VENUS SWINGBY WITH EARTH-TO-MARS TRANSFER ANGLES GREATER 

THAN 180° - Continued 

(b) Date of ideal alinement is November 5, 1981 (Julian date 2 444 913. 65} 

OKBlUL SUrllMt * b.iJL' H'L h l»PSli .LMIuDt - tbt.UU 

INCOMING K-lNf INI Ty* l|V/b*L)|) PA* - l 7 M • / £ L/ h L * J J . (J 2 

OUTGOING V »|Nf 1 N I T V s l] Jil ,uu - ] 7 2 . U b U t. L « |S.JV 


N H LCC |MCL U m t G A Pt«juy Itl uV 1 H OV V p trtl APO AU I 

1 1*U • 7 32 * S 7 4 jb • t> < *Z77* jy.VsH,, 169.7^9 ^6. j |98 V.pjjjd H« u j» 7b2o. 3 2 b 3 j • J6V/6# 


T JO 

NGut 

OMt G A 

1 H U »N 

T M 1 0 V 

b°« 0 V 

V MAG 

ALP A 

OLLTa 

PC mN g 

28889u3 , 6 b 

1 7 J , 2 6 

23 1 .*>« 

7 . b U 

L 6 2 3 . *46 

! 3 1 bo .99 

1 3Md 7 . Uu 

-1/7,7 6 

1 U • 1 3 

2 1 8 9 

Z8889o«.6b 

1 7 i * b 1 

233. 9 f , 

2 . UO 

SUP’ * 03 

1 26 1 7 . 38 

I27uu*uu 

“ 1 7 3 • b 6 

1 1 • 3 i 

1 4 1 S 

28 8*9 | j. 60 

1 6V • 7b 

236 . 1 ? 

■ un 

*4 b n 2 • 6 *4 

1 2330. 1 * 

1 1 V 7 0 • uu 

- 1 / M • 22 

1 3 • J 2 

4 UU 

288*9 | 8 . 6 b 

168*00 

238 . 3u 

• UO 

M 6 „j .22 

1 2 1 30 . 7 J 

1 1 3*4 1 . Uu 

- 1 7 2 • u :> 

l b • 39 

* UM 

298*9^,). 6 S 

1 66 • 2*4 

2 *♦ 0 « b A 

8 . bO 

•4 b *4 6 • 2 H 

l20 7 3.7i 

1 (J 8 1 6 . uu 

“ 1 7 1 • 7 b 

1 8 • bb 

1 4b J 

299*92B . 6 b 

1 6*4 • m9 

242*78 

I 7 • bo 

4 V 7 7 * 9 fJ 

1 2 b 0 b .80 

lUMM 2 *Uu 

- 1 Xu*** 

22 *6 i 

3 * 7 1 

299893;) , 6 b 

1 6 2 • 7 3 

2Mb .On 

3b . 00 

t, H 9 2 • 8 V 

13820.80 

1 U 2 7 6 • uU 

• 1 7U . | 6 

27*68 

6 . 28 

288*938 » 6 b 

1 6 u • 9t> 

2*4 7 . 2 | 

b* . pu 

7 Ub 2 • 02 

1 8 b 7 V , S 3 

1 UMOb * Uu 

- 1 69 . Jb 

33 • «u 

9 t b 3 


* M M A tCC 1NCL * 0 L> t 0 M f G A P t P I 0 U T L 1 U V UK uv V p t*J APO ALl 

2 3 • U 438611*4408 • *t 1 6 8 J 2 V/.U20n ‘♦•1*461 2 2 • S 3 7 J . *, 2 9 *4 7b7^. H 7 5 tt • 29/6J. bbb9. 


T JO 

‘'‘Out 

OMt G A 

1 KU A N 

T M l 0 V 

SUM [jV 

v MAG 

ALP A 

0 1 U T A 

P C ANG 

2 *4 *4 *4 9 Q 3 # 6 b 

2*Ub 

3 U • 3 * 

- 1 . bO 

o 1 m M . 6 b 

1 29U2 . 63 

l JMb 7 • UU 

- 1 / 7 * 76 

1 U • 1 3 

1.6? 

2MMM9(Jb.6b 

3*10 

26 * M i 

- • bO 

7 8 j 8 . q 9 

l 2 b 7 6 • u 7 

1 2 7 UU • Uu 

- 1 7 b . 66 

11*31 

4 60 

2 *4 M *4 9 1 3 *6b 

*4 * 1 b 

22 • Me, 

• on 

7b72* 16 

I 2330 . 1 m 

1 1 V 7 o • uu 

“ 1 7 M , 22 

1 3 • U2 

• UU 

2 M M M 9 l 8 .6b 

b* 1 9 

1 0 . bn 

• un 

7372*71 

1 2 1 3U * 6 V 

1 13*4 1 • UU 

- 1 7 2*6s 

lb* 39 

4 U2 

2MMm 92 j ,6b 

6 • 2 *4 

1 M . b M 

• 1 • UO 

7 2 ? b . 3 3 . 

i 1 963 . 32 

1 U« 1 6 . uU 

- 1 7 1 • 7b 

l 8 *bb 

.63 

2MMM928 .6b 

7*28 

1 0*b9 

“ 1 • bo 

7 2n 2 • 79 

1 1 9 6 U • 7 6 

1 UMM2 • Uu 

• 1 7u , 8 V 

22*61 

1 *VM 

2MMM9 3 j,6b 

H.J3 

6.6m 

• 1 • an 

739b • M2 

1 2 1 b3* Mj 

1 U 2 7 6 • u u 

“ 1 7 u • 1 6 

2 7*60 

3*06 

2*4*4 M9jb. 6b 

V . 3 H 

2 -60 

• 1 • UO 

7/7 3*8(1 

1 2b3 1 *bV 

1 U*4Ub • Ub 

- 1 6 9 • j b 

33 * bj 

8.91 
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TABLE m. - THE 1981 VENUS SWINGBY WITH EARTH- TO- MARS TRANSFER ANGLES GREATER 

THAN 180° - Concluded 

(c) Date of ideal alinement is November 10, 1981 (Julian date 2 444 918. 65) 


0 Kb 1 TAL S 1 At T 1 ML « 

S • (Jit 

rt h l »p si i «Lili nut « 

/b2.UU 





1 N C 0 H 1 N ij ^ - 1 fa Y 
buTfaUINti V.lNf 

IM T y* 
I N 1 T r* 

l | J 4 1 • On 
i 0 a i 6 , i j u 

N<* - | 7 2 >bb ut L* 

K 4 * - | 7 1 , 7 b i, t L - 

l 3 • iV 
l a .bb 





M Sh A 

MO •b49Vdl 4 . U H 

tct 

• b’ub 

iNCL 
i 1 • Jilin 

NllUf 0m III 

168.69/1 2 J 7 . h 71 3 

6 ■ yUUU 

ifcl uV 
b 6 u U • 

C 1 * uv 
6b JU • 

V l't«l 
3|B33. 

4 P 0 *Ll 
1 0*03. 


T JU 

NUlIt 

OMtbA 

T hu A N 

T M 1 DV 

SUM (/ * 

V Mau 

A<■^ A 

ut L I a 

PC A Nti 

244<*VU J» 6 b 

I 77. MV 

226. 7ft 

| 7 . UO 

/ 346 .S3 

1 368 6 • / i 

i J 4 d / • uu 

- 1 //. /6 

i ci • i j 

b • «« 

2HMM9QB . 6 b 

1 /l.UV 

2 3(J • 4 i 

6 • b(J 

64? 9 • 9 9 

1 2 V 6 u • l > 

1 2 7 Uu • UO 

- 1 7S . B 6 

11 • 3 l 

3 t » 3 

24449 i j. 6 b 

1 > 1 > 8 » 

234 . U4 

2 • b (J 

b b 6 2 *04 

\ lilt. it 

1 1 9 / d .uu 

-171.22 

1 3 • lu 

1.39 

2 4 4 4 9 1 B • 6 b 

| 6 b .V[J 

237.6; 

* On 

bbnO • 9V 

12130.6/ 

l l 34 i • UU 

-I72..S 

!&• JV 

• uQ 

244492 3 . 6 b 

1 PS • ?U 

29 1 . 3 , 

• 00 

4443.8b 

1 | ’ 7 ** . U 3 

1 JtJ i 6 . UU 

- 1 / I . 7 S 

1 a • bb 

• 1 U 

24 4 492b . 6 b 

1 Pi • Vo 

2 4 4 • 9 y 

b • 0.0 

b 4 A 6 • 3 1 

1 |VV 6 .t1 

1 U4H2 . UU 

* 1 /u.bV 

22*6 1 

l .92 

244493 3*6b 

|SV.9| 

2 48 . b 7 

1 1 • b(j 

a U 4 3 • 8 2 

| 2 &Z 1 • 0 J 

l U 2 7 6 • u u 

“ 1 / 0*16 

2 7 • 6 a 

4 • d 3 

2444938.6b 

lih.Vl 

2 b 2 -2.1 

3 3 • bfl 

7 1 h 1 • 8 8 

136/2.06 

l U 4 U b . Uu 

“ 1 6 V . 3S> 

3 3 • 8 y 

7 • b D 


M 

SM* 

LCC 

1 NCU 

N 0 U t U rl ► l» A 

Ht.K I UU T L I uV 

Cl* U V V Pl«I 

A P Cl Al-T 

1 • 0 

. H jVl 7U»U» 

.48761 7 

1 2 7 • 7 b Hn 

199.4*3© • u, 8 6 9 

4 • ?8b7 66 Jtt . 

b 4 9 3 • 3U sv b. 

73 16. 


T JU 

N OUt 

0 M 1 (,* 

T HU a N 

T n J Dv 

SUM ijy 

V«Ao 

ALL A 

Ut C I * 

PC ANu 

24449Q3. 6a 

1 d / • 73 

^ 2 5 « 2 1 

1 3 • yu 

7b4U • 4 4 

i JU 3 J . UO 

1 34 a 7 . uu 

-1/7,76 

l U * l 3 

2.61 

24449ob.6b 

1 V j .64 

2 2 8 • 0 o 

7 . uO 

7149.58 

1 i 6 H i . 1 M 

1 2 7 UQ • UU 

- i / b . a 6 

11*31 

1*63 

24449 1 3 • 6 b 

1 Vb.bb 

2 30-79 

2 • bn 

6 b c, o . 7 1 

1 2 3 1 J . i / 

l 1 9 7b . Uu 

- 1 7 4 . 2 2 

1 3 * U2 

.69 

24449 )B. 6 b 

1 99 . 4 6 

233. b9 

• UU 

46*8.1 1 

1 2 1 3 0 • 6 / 

l 1 34 J . UU 

- 1 / 2 . rt b 

1 b • 39 

• UU 

2444923.6b 

203 • 3? 

236 . 3m 

• un 

A 4 A 1 *31 

1 | 9 7 J . as 

1 U* 1 6 . uu 

-1/1.73 

1 b • bb 

• J 8 

244492b . 6 b 

207 • 2H 

239.1 7 

3 • 1)11 

4392 *6 1 

1 1 66b . 1 / 

l U442 *UU 

-1/0.89 

22*61 

. d 1 

244493 J. 6b 

2 11*19 

241*9; 

1 1 • UO 

4 b 4 4 ■ 7 7 

I 203 / . 33 

1 0 2 7 6 • UU 

“ 1 7 U • id 

27 • 6 b 

2 * / 1 

244493b .6b 

2 1 b* 10 

241.7* 

2b .,in 

7 1 AQ • 32 

l 2602 . BV 

1 04ub • uu 

- 169 .Jb 

33 • ©u 

b* 20 
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TABLE IV. - THE 1981 VENUS SWINGBY WITH EARTH-TO-MARS TRANSFER ANGLE LESS 


THAN 180° 

(a) Date of ideal alinement is December 20, 1981 (Julian date 2 444 959. 31) 


ORB I I AL S 1 A 1 T I Mt * 

INCOMING V - INK IN 1 T 1 • 
outgoing v- i N f i n i t r ■ 


S.Ug PfcW I A P S I b AL.TII.lUE » 2*2.00 

10990.00 Ra* l & I « To UtC* 
l a ft 2 9 .00 HA* 199. 3n ufcL* -13.77 


N M iM* t-CC 

1 l.Q .3b2369*U6 .3*U9b9 


1 N C U Nnut 0 m f. G A 

99.6030 l/U.2,0b I«b.69l2 


Hfcfl | 0U I t 1 OV C 1 K UV 
3 • q 7 6 9 7 7 uu . 9 I 6 b . 


V T'EhI 
29160. 


T JO 

NODE 

Out Ga 

2444?H4 . 3 1 

192.12 

163.9* 

2444*49 • -3 1 

184.81 

17 1.19 

2444*54 • 3 1 

1 1 f.b\ 

178.9, 

2444*59.31 

1 >u*2l 

1 8 b • 6 4 

2444 V*4 • 3 1 

1 62-9 1 

1 9 2 . 8 7 

2444969* 3 l 

l5b.*l 

200 • In 

2444974*31 

1 48 • 3 1 

2 0 7 • 3 2 

2444979.3 1 

141 *00 

219. bb 

2444984 • 3 | 

1 3 3 • 7 L) 

22 1 " 7 9 


T KUAN 

TmI Ov 

SUM ov 

19.50 

797 8 . 1 0 

1 2 1 93 .63 

1 U *bO 

7 8 1 b • 8 b 

11981.27 

1 • bn 

7« , 0. 1 7 

, l97 b .b9 

• UO 

7 7 ,, 0 • Ob 

1 16 6 6.98 

• UO 

7 ?b9.62 

1 1 9 2 u • Ob 

2 • UO 

7926.27 

,209 1 . 7 U 

7 • bO 

H 9 b 6 • b 6 

I 2621 .96 

15*00 

95q9 . 02 

,3679.99 

23. bu 

1 1030*52 

16196.99 


V M A b 

aLF a 

DEM* 

1 1 09 V . UU 

l 6J . 9b 

“2 7 • Ub 

1 Uto 7 l • Uu 

lb9.3V 

-25 • 3b 

10^1 4 • U U 

ib3« Vb 

-24 • 1 b 

iu44u.au 

1 b 1 • / U 

• 1 6 • b 6 

10*29 .UU 

14V.J8 

-13*77 

1 1 U9j. UU 

14/ .b/ 

-9*48 

1 1 789 *uu 

146*32 

-b* 9y 

1 2*8 1 »UU 

1 4b . 7 6 

-2*92 

I3739.au 

l 4 b • 8 U 

-*bb 


APO ALT 


9 9 9 8. 


PC ANg 
. 60 
1 • 1? 

1 ■ 66 
«U0 
• 17 
1*31 
9" u2 
7.28 
10.82 


n m sm» tec incc node omega 

2 1*0 .393876*08 .928309 197.7690 ll9.bi.36 |79.6979 


PfcRlOU Ttl UV Cl8 OV V Pfc.Nl 
3.4369 6967. 9878. 298b|. 


APO ALT 
b 6 1 b • 


T JD 

NOUfc 

OtffcbA 

1 HU AN 

T Ml Ov 

2999999.3 l 

96.63 

143.4c* 

v .uu 

7 * n 8 • 03 

2999999.31 

106.73 

1 53 .** 

4.00 

7 2 2 2* 1 * 

2999969.31 

112*69 

1*4.48 

4 • DO 

7 1 ft* *22 

2999969.31 

119.69 

175.00 

• 00 

a9« 7 • 04 

2999969 . 3 1 

126.96 

185.5? 

• 00 

704 1*12 

2999969.31 

1 3 3 * 3b 

196.Q1 

1 *00 

7 2n 1 *85 

2999979-31 

1 9 0 " 26 

20* • 55 

3 • 00 

75ft9 . 99 

2999979.3) 

197.16 

217.07 

5 • 5 n 

«2 1 7 .38 

2999989 . 3 1 

1 69 • 07 

22 7 .59 

* • 5o 

9 0 S 5 .85 


SUM ov 

v mag 

m*»F A 

DE l Ta 

PC AN& 

12986.9/ 

l 1 U 9 9 .00 

163.95 

- 2 7 • Ub 

4.11 

|210U.*U 

1 067 1 .uu 

159.39 

- 2 b • 3b 

2**8 

12069.66 

10*19. UU 

1 b 3 . 9 b 

-24 • 1 b 

2.44 

1 1 6 6 b . 9 8 

10990. OU 

Ibl .70 

— 1 8 * b 6 

• 00 

11919.66 

10629.00 

149.38 

-13*7/ 

• u* 

1 208U .28 

1 1 09o*uu 

147. b7 

-9*48 

1 *03 

! 296b . 93 

1 1 789 .UO 

14*. 32 

-b*9y 

2*89 

13096.81 

1266 1 .00 

1 4b • 7* 

-2*9? 

4.72 

1 3939.28 

13739.00 

1 4b * 80 

-*bb 

6 . *0 


N M SMa 

ECC 

1 NCL 

NOOt 

3 3 ' 0 .369Vbu*08 

. 3 8 2 9 9 3 

8 1 • 1 Buo 

32* .7,3* 


Omega 

b u ■ 9 1 1 6 


PfcN 1 Ou 

3.2932 


Tfcl uV 

7 4 t>b • 


C I « UV 

HHUI • 


V PfcNl 

299 yi . 


APO alT 
9863* 


T JO 

NUOfc 

(j A 

TKUaN 

1 91 Ov 

SUM 0 v 

V Mau 

A^f A 

UE T A 

Pc aNj, 

2999999.31 

333-20 

*3-34 

;Q . 00 

9099.03 

1 34b V . 6 d 

1 UU 9 V « ou 

163.95 

- 2 7 • U 5 

7 .65 

2999999. 3 1 

33 1 . 7 1 

59.01 

b.50 

9129.92 

l 252b. 56 

l U*7 1 .uu 

lbV.jV 

-25*35 

4.9* 

2999969.31 

330 " 2 1 

54.77 

3 • 50 

7620*32 

1 l 9 2 □ • V 7 

1 0* 1 4 • UQ 

1 53 • Vb 

-24* lb 

• 4 U 

2 99 9 969 . 3 1 

328 . 7 1 

5 U • 4 1 

• ua 

7 9 a 9 . 8 3 

, 1 b*b . 4b 

iU44U.UU 

Ibl „ 7o 

- 1 B* 56 

• UU 

2999969.31 

327.22 

4 6 . U, 

• UO 

76,8.83 

1 I'M V. 47 

1 U* 2 V • ou 

|4V.3b 

-13*77 

• U3 

2999969 . 3 l 

32b "72 

4 1 . b n 

2 • b U 

7 6 a 3 . 8 u 

1 2064 .45 

1 l 1 >V u • UU 

147.5 7 

- 9 • 4 d 

• 6 3 

2999979-31 

329 .22 

37.49 

7 • 5u 

7961 .6U 

l 2352 • 24 

1 l7eV.au 

1 46 . 32 

- b * 9 a 

1 . b 5 

2999979.31 

322-/3 

3 3 . 1 M 

l * *Q0 

99*2.37 

l 2«*3. 02 

12 * 6 1 • Uu 

145.76 

-2*92 

3.5b 

2999989 . 3 1 

32l .23 

28 . b ; 

25 • bn 

9272*69 

1 3*23. 34 

1 37 3V • Uu 

1 4b • BU 

-*bb 

5*54 
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TABLE IV. - THE 1981 VENUS SWENGBY WITH EARTH-TO-MARS TRANSFER ANGLE LESS 

THAN 180° - Continued 

(b) Date of ideal alinement is December 25, 1981 (Julian date 2 444 964. 31) 


OSBlTAL STaTTIME * 

5 • Ou 

PE h I»PSl S alTI lUDt » 

262.00 



INCOMING V-iNF INI Tf» 
OUTGOING v.ineinitt. 

l 0t>29 «un 
1 1 U90 >on 

RA» |49.36 OtL* 
SA* |47,S7 utL* 

-13.77 

-9.48 



M SMa ECC 

• 0 «37|677 + OA • 39 4 j 6 | 

1NC.L 
93* 1500 

node Omega 

164.5348 I 6 1 . 7 b 2 b 

per i oo 
3.3333 

Tel uv 
7900. 

CIS OV V PfcH 1 APO ALT 
9b 1 9 . 29^22 • S0B4. 


T JO 

NODE 

OMfcG* 

tnuan 

1 HI DV 

SUM ov 

v 

A ^F A 

de l ta 

PC A N b 

2999999 . 3 | 

|9|.«2 


30.00 

79a 3# 8 1 

12982.93 

11099.00 

163.95 

- 27 • OS 

3.89 

2999999.31 

189.00 

1 68*9a 

I’.bO 

7789.76 

r23o8.9o 

10671 .OU 

159.39 

"25 • 3b 

3.78 

29999b9.3l 

1 7*. 1 7 

176.22 

8 .SQ 

7 7 aH • 58 

12283.71 

i06iH.au 

153.95 

•29 • l S 

3.9J 

2999909 . 3 1 

1 7 | * 3b 

iaa.9«> 

2.5n 

7397.31 

1 1 9 16.93 

IQHHO.UO 

1 b 1 . 7 0 

• l 8 • b6 

1 .98 

2999969*3 1 

I 69 «S3 

191.7c, 

•00 

79 0 q.32 

1 19 19.99 

10629.00 

l#9. 38 

-13*77 

.00 

2999969.3) 

1 b 7 • 7 1 

199.52 

• 00 

7536. 13 

I20SS.26 

i 1090*00 

197.87 

-9*98 

• 19 

2999979.31 

| bo . 89 

2q7 . 2* 

2*00 

7010.92 

12357. b9 

1 1 789.00 

196.32 

•b • 90 

1.99 

2999979.31 

1 99 • 07 

2 1 5. US 

7 *50 

8 5 14 8 .95 

1 306B.Q7 

1268 1 .00 

195*76 

-2*92 

9.78 

2999969.31 

137.25 

222*8! 

15.00 

9 7 H 6 • 8 a 

1926b. 93 

13739.00 

1 9 5 . 8 U 

-*bb 

7.9Q 


N M SM* ECC 1NCL 

2 1*0 .451643 L&u»|86n 


node 

124. 0*02 l rt 3 • c. l 6 3 


PtRlOU Ttl UV CIS UV 
3.871U 6796. b 1 2 2 ■ 


V Peal APO ALT 
3Q 1 24 • 6367 . 


T JO 

NOUt 

0«L&» 

tkuan 

TM i DY 

SUM OV 

VMAG 

alf a 

dei-Ta 

PC aN6 

2444944.31 

90*29 

192*97 

5- 00 

A9s9 *99 

I 208 I .02 

1 1099.00 

163. 9 * 

"2 7 • 05. 

• 87 

2999999 . 3 1 

109.73 

15 2.71 

2 *00 

a8 | 7 • 73 

1 1939. 3u 

^067 1 .UO 

159.39 

"25*35 

• 52 

29999S4.3I 

1 1 1 • 1 e 

162*99 

8.50 

7277.27 

1239U.8S 

‘l Q614.0U 

153.95 

"29» 15 

9*23 

29999SV.3I 

1 17.62 

173*26 

1 -00 

6 7 a 6 . Q6 

l 1 9 q 7 .64 

10440. UU 

1 5 1 • 7u 

"18.56 

1*20 

2999969.31 

1 2H • 06 

183*5? 

•00 

A 7 9 7 * 8 6 

1 ) 9 ) 9 » 4 4 

10629.00 

199.38 

- 1 3*77 

• 00 

2444969 . 31 

1 3q.50 

193.7a 

•00 

A9 1 3 • 5 2 

12055.09 

1 1090.00 

197.57 

•*9*98 

• 10 

2444974.31 

1 34 .9h 

2 0 9 • 0 h 

• 5n 

71 A5.87 

1 2287.4b 

l 1789.00 

196.32 

-5» Vy 

• 85 

2444979.3) 

193-39 

2)9. 3n 

• 5 0 

7Sc.3 -59 

12675.12 

12681.00 

195.74 

-2*92 

2* 10 

2444984 . 3 ] 

199.83 

229.5a 

•50 

ft 1 33 .3 1 

1 3254.89 

13739.00 

195.80 

-•55 

3 • 7 1 


M SM* 


ECC 

1 NCL 

node 

Oh . ga Period iei uv 

C 1 * 0 V 

V PfcKl 

apo act 

J *0 • “ l ’ ' 1 9 + 00 

H 6 35t)f) 

8 1 • 0 7 S n 

327.1743 

46 . (16 1 9 

4.UUU0 667b. 

5299 . 

30*97. 

6665. 

T JO 

n o a E 

0 *EGa 

TKUAN 

T Hi DV 

SUM DV 

v Mag 

a l f A 

D£ l Ta 

PC aNg 

2444944.31 

331.56 

5 8 . 5 n 

29.00 

4 6 7 3 * 5 U 

13917.86 

1 1099.00 

163. 9b 

-27*05 

8.9 1 

2444949.31 

330-96 

55*39 

1 7 m UO 

7 6 n 6 • 7 2 

1 26S 1 .07 

1067 1 .Qu 

159.39 

-25*38 

5.7b 

2444954 .31 

329*37 

52.20 

9 • Sn 

A7 n 4.92 

1 1949.27 

10614.00 

1 53.95 

"2 9*j5 

• 8 1 

2444959.31 

328 . 27 

99.17 

3 • 00 

66 3 b . 02 

1)879.38 

10440*00 

l 5 i • 7 u 

-18*56 

.69 

2444964.31 

3 2 7 - 1 7 

96 * 0 A 

• 00 

*6 7b . 09 

1 1 9 1 9 . 44 

10629.00 

199*38 

-13*7 7 

• UO 

2444969 . 3 1 

326*00 

92 • ’5 

• 00 

A 8 | 0.« 4 4 

1 2054.60 

1 1090.00 

1*7.57 

-9 *98 

• U2 

2444974.31 

J2h *96 

39 .8h 

3 - On 

7 0 1 8 . 2 7 

1 2282.63 

11789.00 

196.32 

rS • 9 (j 

• 67 

2444979.31 

323*88 

3 6 • 7 i 

’*00 

7432. S2 

1 2676.87 

1 266 J .00 

195.76 

-2*92 

2*07 

2444984.31 

322*79 

33*62 

1 T .50 

6056 . 76 

1 3301 . 14 

13739.00 

1 9b • ttj 

-•55 

3.77 
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TABLE IV. - THE 1981 VENUS SWINGBY WITH EARTH-TO-MARS TRANSFER ANGLE LESS 

THAN 180° - Concluded 

(c) Date of ideal alinement is December 30, 1981 (Julian date 2 444 969. 31) 


o n bual sTArrlMt • 

s .ol. 

PE k 1 a P5 1 S altITijDE » 

262 .GO 



INCOMING V-iNFIMtt. 
00TG0IN& V-InF IM I r > 

1 1 090 . OU 
j I 709. oo 

r A« |97.57 DfL- 
«A« 1*44.32 litL- 

• 9 • 40 

-b • 9|J 



M SM A tec 

I • Q .910*39.08 .95 1693 

l^CL 
3 * *026n 

NCOb 0h t (,A 

149.9*2’ 1 4 8 • i 5 3 1 

PE« l Ob 
3 • b / 1 U 

It 1 U* 
6 9 J 3 • 

tIK DV V Pt-Kl APO ALT 

b 1 22 . 3Q i 2 4 • 6367* 


T JO 

NODE 

0M t 

T kuan 

T mI 0 V 

sum D v 

v"ag 

A L F A 

UE L TA 

PC ANb 

2*4*4*4.44*4. 3 1 

1«8.30 

1 6 U . 7 n 

4b • 50 

0b (0 . 27 

13431.8b 

1 1099.00 

103. 9b 

• 27 • Q b 

7.97 

2 * 4 * 414 . *49*31 

| 6 2 • & *4 

160.19 

34. oo 

02 74 . 4 l 

1339S. 99 

10471.00 

1 S 9 • 3 9 

“25 • 3b 

7 • bi 

2999*59 . 3 1 

174.77 

1 7&.6b 

2 1.50 

0160*32 

13287. 9u 

1 04 1 *4 .00 

1 b 3 • 9b 

“29 • 1 b 

MS 

2999*59 * 3 1 

1 1 1 *00 

103.17 

1 1 *00 

7277*06 

[ 2399. *4 *4 

10*4*4 0*00 

161 • 7u 

" I0*b6 

9 .90 

2*4 *4 44.4*4 . 3 1 

1 4b . 23 

l 9 q . 66 

2 * 50 

6003 •‘♦8 

I200S.QS 

10*29.00 

199.30 

“13*7/ 

1 • 0 3 

2*4*4*4.49.31 

l S 9 • *4 4 

1 90 . is 

•OU 

6 9 j 3 • 2 0 

12059.78 

1 1 09u • 00 

197* b7 

-9*96 

• UO 

2*4 *4**. 7*4 . 3 1 

1 b3 . 49 

2o & ■ 6 h 

•00 

7 i 14 B • 0 3 

I227U.01 

1 1 749.00 

196*32 

-5»9U 

• 06 

2*4*4*4.79.31 

1*47.9 3 

2 1 3 • 1 1 

2 • QO 

74*9 • &9 

12*11.17 

1 244 1 . 00 

19b* 76 

-2*92 

1.95 

2*4*4*498*4 . 3 1 

197*1* 

220^7 

7 *U0 

8 | SH • S 7 

1 3274. IS 

13739.00 

1 9& • 00 

“•bb 

3.03 


M SMa 

1*0 •*429 3ot>.ua . 

ECC 

9 7 b 9 8 9 

INCL 

1 b 2 • 1 8 S 0 

node 

12 ’. 1 1 1 7 

OMt&A PERIOD Tfcll 

[92.49S2 9.(379 6*8 . 

.V tl R DV 

. 5344. 

V 3 0i^‘. 

apo alt 

4*8 | . 

T JD 

NODE 

OMEGa 

T K U A N 

T M I OV 

SUM OV 

V "AG 

A l Fa 

DE L Ta 

Pc aNg 

2*4*4.9*4*4.31 

99.38 

193.71 

7 .on 

* 898 .40 

12224.40 

1 1099. OU 

143.95 

-27-05 

2.99 

2*4*4*4999.3 1 

10S.33 

ISj.Sn 

* « 5n 

A*n 2 • S l 

1 2270.51 

10*71.00 

159.3* 

-25.35 

3.97 

2. ..959.31 

111-27 

143.2. 

1 3 .50 

7S73.S2 

I299 t ,S| 

10*19.00 

153.95 

-29-15 

5.98 

2999*59 .31 

117.22 

173-07 

4*00 

*7*1.19 

12129.18 

10990.00 

151.70 

-18.54 

3 • u* 

29999*9.31 

123-17 

182 -8* 

• 50 

A b u 3 * 90 

11*51 .90 

10*29.00 

199.JB 

-13.77 

l " 0 * 

2999949.31 

12*«11 

192.4s 

•00 

*4.4.78 

12059.78 

1 10 * 0.00 

197. a7 

-9 • 40 

.00 

29.9979.31 

1 3b. 04 

202 -. » 

• 00 

* 9fl 1 « 9 4 

1 2249.9* 

1 1 789 . OU 

1*44.32 

“5 • 9(j 

• 03 

299. 979 . 3 1 

191.00 

212.22 

“ 1 »00 

72n2.89 

I 2 5 7 0*84 

12481.00 

195.7* 

-2*9* 

.*1 

29.9989 • 3 1 

1 9* .9b 

222.01 

“2 • 00 

7*27.04 

1 2*95.0* 

13739.00 

195.40 

“•bb 

1-81 


M SM A 

tec 

JNCL 

node 

0 N p Ci A 

PtKiOu T E 1 uV 

ClK OV V P E k 1 APO ALT 

3*0 *5323004-00 

• 5 7 7 0 3 7 

0u • 9 1 9 q 

3 i 0 • 0 M g 1 

4 2 • ©25 7 

5.7143 b 6 b 9 • 

*394, 3 1 39d . 10379. 


T JD 

NODE 

OMt G a 

TRUA N 

T Hi DV 

SUM ov 

V MaCj 

A L F A 

delta 

PC 0N(, 

2999*99. 31 

329 • «*u 

52.2s 

4 6 • 0 n 

7930.66 

1 9326 . 37 

1 1 U99 . UU 

163.9b 

•2 7 • u 5 

9*10 

2999*99.31 

320* 73 

50 - 3*> 

36 • on 

6065*09 

13260.81 

1 U67 1 . Uu 

159.39 

“2b«3b 

6*40 

2999*59.31 

320*06 

40 . b2 

24. UO 

b 7 s 6 • Jo 

l215l.8i 

1 U6 1 4 • UU 

153.9s 

“24 • 1 b 

2*16 

2999*59.31 

327.39 

46.66 

l 2 • bO 

5 b 9 6 * 2b 

1 1*9 l .*/ 

10440. QU 

lb 1 *7U 

“ 10 *56 

1.79 

2999**9.31 

320* 7 1 

44*79 

4 • U 0 

bb 4 2 • 4 4 

1 1*34.16 

l U6 2 9 • u u 

14V.J0 

“13*77 

. 7 1 

2999**9.31 

326 *04 

42*93 

•OU 

b 6 S 9 • o 7 

1 205 9 • 7 8 

1 lOVj.Ub 

147*57 

-9 . sd 

• UU 

2999*79.31 

32b- 37 

« I -06 

• 0*' m 

5074*23 

1224* .9-4 

i l 7 6 9 • 0 -j 

1 06-32 

-b*9u 

* u 1 

29*44*79.31 

320*69 

39 . 1 0 

3 • b c 

A 1 «2 • 9 4 

12588.46 

1 260 J . Ob 

1 Hb . 76 

• 2 • 9 2 

.92 

2999*89 . 3 1 

324 *U2 

37.31 

J U *UQ 

6 6 9 9 *Q6 

1 3 0* 9 . 7 8 

137 3 9 . ij u 

1 4b . bu 

-•bb 

2* 30 
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TABLE V. - SOLUTIONS EVALUATED FOR THE TWO-IMPULSE DIRECT MARS STOPOVER IN 1986 

(a) Type I heliocentric transfers 


ORBITAL STAfTlHE » 

& .C0 

PEKIAP51S *LT!TUOE 

■ 2 & 2.00 



INCOMING V*INF|NITY« 

out go i n g 7 v - infinity* 

%7S *00 
9537.00 ' 

«A* "63*59 u EL ■ 

RA» -67. 29 DLL* 

“52*78 

-49.47 



N -i $rt A ECC 

INCl 

• node 9*jeg_*_ 

PENIwU 

iei pv 

cr*_yv. v pfcH j_ a p_o a i , t 

I 1*0 .2532*9+08 • ! I Q9 \ 7 

55 • 6 | 6 c 

.6903 " 134.2975 

| . 0 75o 

10259. 

135q# 26353. 1107. 


T JD 

4006 

omega 

THUAN 

.. .Tin J?y . 

..syn_D9 

¥M_AG 

ALFA 

. .oecja 

. -P.c ang. 

2996520.00 

79.60 

95.39 

28.50 

1 ?763. 46 

19113.61 

1 1 225.00 

-42.35 

-6 1.69 

1 1 .7o 

2996521.00 

59.82 

103-12 

26.0J 

12995. 19 

13895.29 

10950.00 

-93.31 

-62.92 

1 1 .09 

2446530,00 

95.03 

110.90 

17.00 

1 1557.90 

12908.06 

10336.0c 

-43.69 

-60.92 

0.34 

2996535.00 

30-25 

118.68 

’•00 

10& 13. 19 

12163.33 

9007 . 00 

-5+.0S 

-50.56 

S.9| 

24H6590.00 

15.47 

126.47 

2.5C 

1 p367 .56 

1 1717.71 

9600*00 

-59. 10 

-55.05 

2.36 

2996595.00 

.69 

139.25 

• QQ _ 

_ln2Sli77_.. 

. --UAJIO.’Z 

9925*00- - 

-63*58 

-52. 7#.. 

»ou 

2996550.00 

345.9 | 

14203 

• 3Q 

10262.07 

11*12.22 

9507.00 

-67.28 

-49.47 

.02 

2996555.00 

331.13 

199.9! 

3.00 

1 0 5 2 1 .92 

1 1672.07 

969 | .00 

-70.12 

-46.46 

3*78 

2996560.00 

316.35 

157.59 

1 1 .So 

1 f 560 .85 

12911,00 

10017.00 

-72. l 1 

-92.73 

e.63 

24*46565.00 

30 1.57 

165-34 

22.50 

13488.04 

15238. 19 

10475.00 

-73.30 

-39.57 

15. OC 

2446570.00 

286.78 

17 3*16 

33.50 

17214.40 

18569.69 

1 1 050.00 

-73.79 

-36.66 

22*35 

2446575.00 

272.00 

1 81. 99 

9 9.00 

J L226 *A3 

. 225 7 6.. 89 

1 1/3J..QO 

- 73. .67. 

09.318 

30.87 

2 *46500. 00 

2S7.22 

189.72 

57.00 

26 1 5 7 . 23 

27507 . 30 

12259.00 

-70.93 

-31*12 

92.1* 

2996585.00 

242*44 

196.50 

77.5o 

31338.97 

32^09. 12 

12552.00 

-65 .05 

-27.U0 

56.5c 

2996590.00 

227.66 

209.28 

1 1 O.So 

3995 3 .36 

36300.52 

12862.00 

-60.88 

-22.3 1 

7Q. 24 


N M 


ECC 

INCL 

j^o p f 

0 _MEG\_ 

periuo 

Tt I DV 

Cl* O’. 

V PER L 

. AP.fi ALT. 

2 3.Q 

.378729*08 

.905933 

77.933 U 

1 00.0745 

83.54 15 

3*4206 

6966. 

.638. 

29641 . 

531*. 


. -LIB 

NOPE 

ONE&A 

J iiiLAli 

LH4 -DV 

-SUM _0¥ - 

¥Ma5_ 

-AMT-A- 

liEUA 


2996520.00 

1 J9.42 

101.01 

.00 

7749.94 

12388. 17 

1 I 225.00 

-42.35 

-6 1.69 

3*46 

.29 9 6525^0 

. ID?.55_ 

97.52 

6 . 5 0 

7576. 13 

12219.35 

lO’bO.OO 

-H3*3» 

- 62 *?^ 

2*09 

2996530. QO 

105.60 

94.02 

4.5c 

7320.03 

1 1966.25 

10336.00 

-49.69 

-60.92 

2*31 

2398535.00 _ 

loa.ar 

90*53 

2 . So 

7120.33 

1 1746.56 

9887.00 

-5^.05 

-50.50 

J .49 

2996590.00 

10 1 .94 

87.04 

.50 

7007.06 

l 1695.29 

9600*00 

-59 • I 0 

- 55 * 0 b 

• 6 1 

2996595.00 

IQnt07 

— 0-3.5^— 

fcUO- - 

^9^.t_^6-9. - 

1 1*03.9-2 

997-5.0Q 


-*2-7* 

w-eti 

2996550.00 

98.20 

80.05 

*00 

6 9 7 4 • 0 7 

11612.30 

9507 .On 

-67.28 

-49.47 

• 05 

2996555. nn 

95*3 9 

76.55 

2.00 

7034.35 

1 1672.58 

9*9 1 .00 

-70.12 

-46.46 

.70 

2996560.00 

99.97 

7 3.-06 

6.00 

7267.73 

l 1905.96 

10017.00 

-72.. 1 l 

-42.73 

2.51 

2.194565*3)0 

92.40 

6 9.56 

12 . So 

7731.04 

i?369.27 

10975.00 

-73.30 

-39.57 

4.45 

2996570.00 

90.73 

66.07 

2 1 . GO 

0434.21 

13072.99 

1 1050.00 

-73*79 

-36.60 

6.52 

2996575. tip 

8 8.36-. 

A?*^? 

_ - 

— 933*^96- 

- 13979.31 

1 1733.00 

-7S.67 

-J_4*O0- 

8 -» 7 3 - 

2996580.00 

06.99 

59. 

43.00 

1 1 05 L • 0 1 

15688.23 

12259.00 

-70 *93 

-31 • 12 

12*99 

2996585.00 

>*•12 

.5 5 1 5 9 

5 7.QO 

1 >639.78 

18278.01 

12552. 00 

-65 *05 

-27.00 

19.50 

2996590.00 

83.25 

52.09 

70.00 

16304.65 

20992,87 

12862*00 

- 60*^6 

-22*31 

26*00 
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TABLE V. - SOLUTIONS EVALUATED FOR THE TWO-IMPULSE DIRECT MARS STOPOVER IN 1986 - Concluded 


(b) Type n heliocentric transfers 


UKB t r*L SU H | » 

b - DM 

p F. w I A P S! S 

At T l 

TUr> p — 

2 6 2 • > J 0 




1 .JC0II1NG V-l ► f 1 lit T»t 
OUTGOING </ - 1 IjF ] u I 1 Y » 

9 9 <* 1 1 . 9*) 
969ft . S » 

9 A = - 9 ? - 

A * -HI# 

* 2 
7 £ 

OFL ■ 
OF L * 

I • 3 1 
-.?l> 




G 1" A t 1. C 

Li . hll 1 nS? + ll<> . bSllbP? 

IMCL 

7/.7R.10 

* • 0 0 F 

1 | 7.r>B7<5 

n M K A 
2 1 1 • 6 1 36 

ptr i ul> 
b-2l 7H 

T t I D V 
5b*6, 

C I 9 n V 
6131# 

V PfRI APO alt 
3113*1. ?3H3. 


T JD 

r.out 

n i-i f (i a 

TrllAlj 

T ►* l Pi 

S E J M P V 

VM*f, 

ALFA 

oel t a 

pc ang 

2HH6H8S. 

31 1 • ii 7 

2 3 3.30 

ii.no 

) 1 HjO-pb 

17962-19 

13011 -5P 

-Gl .27 

1 *63 

L 3 . 65 

2HH6HV J . j 6 

311.72 

2 3 1 .19 

l n . b u 

10695.9) 

1 6 8 2 7.25 

12670*61 

-G8.8S 

2*50 

11*91 

2446495. 36 

3 i 2 # 3 / 

2 2 9.fl0 

9 . 8o 

9650 • ft ? 

1 b 7 fl J .97 

1 2 3 3 7 • 0 H 

-56 .53 

3.22 

10*29 

2HH6buu# 3b 

3 1 3.1.3 

2 2 6 . 8 a 

»< • bn 

h 7n3 .b*> 

1 H834 . 93 

1 20 | 0 .n* 

-5*4.3? 

3-70 

8.71 

244650b . 3b 

3 1 3 .<>8 

22H .fa i 

7 . ny 

7 « 5 2 . 1 ? 

l39fl3.H7 

1 1 6b_7 • 4 3 

-S2.1S 

4*02 

2*J4 

244ftbl 0* 36 

3 1 H • 

2 2 2 . M j 

4 .SO 

7114.97 

1 3 ? M 6 . 31 

1 \ 37n.3fl 

-50. OP 

4.05 

5.67 

244obl b.3b 

3 1 H * *9 V 

220 .23 

2.10 

6S0S * A ft 

1 2637 .02 

1 1059. n 4 

-**a. l? 

3 • 82 

4. 14 

2446520*36 

3 1 S • 6 H 

2 10.03 

-1.00 

6047.42 

1 2 1 7P. 7ft 

10784.71 

-SA . ZS 

3.29 

2 • b 8 

2HH6b2b.ua 

3 i 6 . 2 b 

2 1 A . n u 

3b7.no 

b 78b. 1 9 

1 1916*53 

10480.66 

-88.45 

2.5 5 

1 *Q7 

2446530*00 

31 6. 9U 

2 1 3 • 6 i 

- 1 .50 

5 6 7 1 . b 4 

11 802 . Hfl 

1 0 2 0 2 • b H 

-83.5A 

2-17 

*34 

244653b* UU 

317-56 

2 11.61 

• on 

5b9S .9t, 

1 l 727.2ft 

9940*94 

— 8 ? . A ? 

1 - u 

• QQ 

2HH6SHIl.au 

3 1 fl • 2 1 

209 • H2 

• no 

b b 3 1 .PA 

1 1 662-41 

9 6 9 8.51 

-8 | . 7A 

-•20 

• 00 

2HH65HS.00 

318.87 

207.23 

-2. bo 

bH 76 .8n 

1 1 607. ft a 

9H 7 9. 56 

— 8 Q » 8 7 

-2*70 

*22 

2 H H 6 550 • UU 

3 1 V • '_> 2 

7. U b • ' J H 

3b 1 .bo 

b 4 4 6 • 4 b 

11577.79 

9298.74 

-39.71 

— 6.86 

• 42 

2446555* OU 

3 2U • | 7 

202 • 84 

33fl.no 

bb 1 3. 1 « 

1 l A 4 H « b 2 

9 2Q 3 • ? 7 

-37. AS 

-14*27 

• 66 


n 5 m#. 


tc C 

I * Cl. 

>(01*1 E’ 

MMF5/J PFRIOd Tf r D V 

CIR nV 

V PpRT 

*P0 ALT 

3*0 .873u41+nB 

742078 

33.0620 

139.335a 

? 7 . 9 3 1 9 12 

•0000 37?9. 

7998 . 

J30Q0- 

21SBJ* 

T JO 

POlifc 

0 M t G A 

T K r A N 

Tm! r.n- 

su m 

VM a6 

AI.FA 

OFL T A 

PC ANG 

24H64Hb* 36 

153.0b 

7.7 3 

59. no 

1 0002-70 

1 8799.99 

1 3 0 1 1 . 5 A 

-61.27 

I * 63 

17*01 

244649U. 3to 

1 M . ft / 

9 . 76 

53. fin 

9/31 • 5 P 

1 7779. 3y 

1 ?6 7 n. 6 1 

-58.8b 

2.50 

14,27 

24464 V{, • 36 

1 bU • 28 

1 1 .Dn 

<•6.50 

8 7 Of' . I A 

1 6 6 9 7 .^6 

12337.0* 

-56*53 

2*22 

11.88 

2446500 • 36 

1 4 fa . 90 

1 3.83 

4u.no 

77 1 3.35 

15711.15 

1 2 0 1 n • n 6 

-54.32 

3.78 

9*78 

2H46SUS . 36 

1 47.52 

15.87 

3 3.50 

6798 . 35 

1 4796. 1 5 

1 1637*43 

-5?. 1 5 

4*02 

7*97 

244651 0*36 

146.14 

17.90 

27,50 

5986.69 

1 3954.49 

1 1370*38 

-50.09 

4*05 

6*32 

2446b l b* 3b 

144.76 

19.94 

2 l . P 0 

5216.2? 

13214.02 

1 1 089.04 

-48 . 1 2 

3 *Q2 

4*87 

2446520. 36 

1 43 . J8 

21.97 

1 s.^o 

4 6 1 0 . 1 ? 

1 7 60 7 • 9 | 

1 r)7 S q • f I 

-46.25 

3*29 

3.54 

244652b . UU 

142*10 

23. Ho 

11.00 

4 1 9 3 . J 7 

! 2 1 9p . 9 ft 

1 11480.66 

-44.65 

2*55 

2*43 

2446530 -UU 

140.72 
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TABLE VI. - THE 1981 VENUS SWINGBY WITH 60-DAY MARS STOPOVER, THREE-IMPULSE 
TRANS- MARS TRAJECTORY WITH MINIMUM-SUM AV (2AV) 

(Date of Ideal alinement Is December 1, 1981 (Julian date 2 444 940. 00) 
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1 NCOM I N(j 
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TABLE VII. - SUMMARY OF EARTH DEPARTURE PARKING ORBITS FOUND FOR 


THE THREE-IMPULSE VENUS SWINGBY MISSION IN 1981 


Epoch, 

month/day/yr 

Inclination, deg 

Apoapsis altitude, 
n. mi. 

Tracking time 
below 12 600 fps, 
days 

11/01/81 

36. 81 
99.02 

23 676 
6 084 

15 

25 

20 

25 

11/06/81 

38.67 

96.84 

13 857 
5 316 

11/11/81 

39.21 

8 050 

15 


130. 00 

5 815 

20 

11/16/81 

39.36 

4 651 

20 


132.61 

3 718 

15 

11/21/81 

40. 51 

2 186 

15 


132.71 

1 752 

15 


96. 10 

2 811 

25 

11/26/81 

77.87 

779 

30 


57.01 

856 

10 


141.20 

3 100 

10 

12/06/81 

75.79 

1 752 

35 


53. 17 

1 752 

25 

12/16/81 

74.46 

2 424 

40 


50.25 

2 074 

30 


159. 10 

3 889 

30 

12/21/81 

71.00 

3 552 

45 


49. 10 

2 811 

30 


161. 50 

4 250 

25 
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APPENDIX 

SOLUTION FOR THE INCLINATION AT THE MINIMUM-ENERGY POINT (v w = 0) 

DURING THE LAUNCH OPPORTUNITY 


Normally, during a launch opportunity, the magnitude of the hyperbolic excess 
velocity goes through a minimum. The solution for the inclination at this time is ana- 
lytic and provides a useful insight into the number of solutions available. It also pro- 
vides convenient initial guesses for the inclination when the complete solution is 
required at some other time during the launch opportunity. 

With the condition = 0, equation (16) becomes 


°oo +f6 oo 


= h 




(Al) 


This equation reduces to 


£ 3 + a 2 i 2 + + q q = 0 


(A2) 


This is a cubic equation in 4 where 


r * 2. 

4 = tan l 


(A3) 


°2 ~ 


8 + tan^^ + ( 2 + s ec 2 6 oo j 

x 


a, = 16 + 8 tan^C^ + —%(% + sec“6„J (2 sec“6 - 1 J 


,) (2 sec 2 6 - l) 


*0 " 


16 tan 1 ^ +^(l - 2 sec 2 6 oo ) 
x 


= a Jh 

00/ o 


(A4) 
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By inspection of these coefficients, a ^ and are always negative and a ^ is always 

positive. This indicates that there are three real positive roots for £. The three 
solutions for | may be written (consult any mathematical handbook) 


where 


h = 2 ^r cos l 


4 2 = 2 ^c°s'(| + 120 j 
h = 2 ^T c°s (f + 240°] 


(A5) 


zil «: 

2 V" 27 


v- 1/2 


< p = cos 


B = i( 2 «2 3 - 9 o'ia 2 + 27 % 


(A6) 


Each part of equation (A5) actually represents two solutions for the tangent of the 
inclination because 



For the posigrade orbital inclination, the solution is 



For the retrograde orbital inclination, the solution is 



(A7) 
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For each solution for i. which is posigrade, there is a solution which is retrograde. 
However, from the restriction given by equation (18), only one is valid; that is, if 
is positive, i can only be retrograde. Thus, there are at most three values 
of the inclination that will satisfy equation (A2). 


NASA-Langley, 1972 30 S— 289 
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